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ABSTRACT 

We put stringent constraints for the first time on the dust properties in the circumburst medium of a 
gamma-ray burst (GRB) afterglow. This is based on the optical spectrum of GRB 020813 {z — 1.255), 
obtained with Keck I LRIS 4.65 h after the burst. From the absorption lines in the spectrum, we 
derive very high column densities for six heavy elements with different refractory properties. The 
relative abundances resemble the dust depletion patterns in the Milky Way, from which we infer a 
visual extinction of Ay — 0.4 and Ay > 0.3 at 95% confidence level. However, the high columns 
of metals and dust contrast with an observed UV continuum spectrum that is remarkably close to a 
power law of the form F^, oc i^""'®, with no sign of curvature, or a 2200 A extinction feature, suggesting 
low reddening. The Milky Way or Magellanic Cloud reddenings are possible only for very low visual 
extinctions {Ay < 0.08 or Ay < 0.2, respectively at 95% confidence), inconsistent with the Ay values 
inferred from the depletion analysis. If we assume a GRB intrinsic spectrum and an extinction law of 
the forms F{ = Fy(5500/A)" and Ax = Ay (5500/ A)"*", we obtain (for Ay = 0.4) the constraints from 
continuum spectrum: 7 < 0.85 and a < 1.72. 

Subject headings: cosmology: observations - gamma rays: bursts - galaxies: abundances - ISM: dust, 
extinction 



1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) are probably 
associated with core-collapse supernovae (Hjorth et al. 
2003a) and can be used to probe the interstellar medium 
of their star-forming environments. One of the open is- 
sues on the nature of GRBs is the dust properties in 
the circumburst medium. On the one hand, high heavy- 
element column densities have been measured in a few 
GRB afterglow spectra (Savaglio et al. 2003; Vreeswijk 
et al. 2004), suggesting that the dust content might be 
important, especially if star formation is very active. On 
the other hand, the spectral energy distribution (SED) 
of other GRB afterglows indicates that dust reddening 
is surprisingly low (Galama & Wijers 2001; Hjorth et al. 
2003b). Up to now, a simultaneous analysis of heavy- 
element abundances to measure the dust depletion, and 
of a suitable UV/optical spectrum to estimate the dust 
extinction law, in the same GRB, has never been done. 
Afterglows fade so rapidly that obtaining high-quality 
spectra, for a full investigation of the dust properties, is 
very difficult. 

The spectrum of GRB 020813, obtained recently by 
Barth et al. (2003), offers us this opportunity. The opti- 
cal spectrum has a very high signal-to-noise ratio (S/N) , 
is very well flux-calibrated, and shows very little devi- 
ation from a perfect power law over a large wavelength 
interval {F^ ex j/-n-92 in AA = 1430 - 4020 A rest frame; 
Fig. 1). Moreover, the absorption lines of numerous dif- 
ferent elements, with different refractory properties, have 
been detected with exceptional accuracy, thanks to the 
excellent spectral resolution. By visually inspecting the 

^ On leave of absence from INAF, Osservatorio Astronomico di 
Roma, Italy 



GRB spectrum, it can already be noticed that there is no 
sign of the 2200 A extinction feature {v « 1.4 x 10"^^ Hz), 
or curvature on even larger wavelength intervals, suggest- 
ing either an extinction law different from the Galactic 
type or a low dust content. 

Our approach to studying the dust depletion is simi- 
lar to what has been done routinely for years for strong 
absorption systems detected in QSO spectra, the so- 
called damped Lya (DLA) systems. The spectrum of 
GRB 020813 shows similar absorption features, most of 
which are associated with low-ionization species, so we 
consider this absorption system a GRB-DLA. Barth et 
al. (2003) detected seven low-ionization species (Fe II, 
Si II, Zn II, Cr II, Mn II, Mg II and Al II), and the hues 
are so strong that the absorbing gas is very likely neutral 
(H « H I). QSO-DLAs provide detailed information on 
the neutral interstellar medium (ISM) of high redshift 
galaxies (Prochaska 2004 and references therein). 

Our paper is structured as follows: in §2, we describe 
the measurement of the heavy-element column densities 
in the GRB 020813 circumburst medium; in §3 we derive 
the heavy-element depletion pattern; in §4 the visual ex- 
tinction is estimated; in §5 we constrain, given the dust 
depletion, the extinction curve; and in §6 we summarize 
the results. 

2. MEASUREMENTS OF HEAVY-ELEMENT COLUMN DENSITIES 

For the analysis of the GRB 020813 afterglow, we use 
the spectrum obtained with the LRIS spectrograph on 
the Keck I telescope in 2002 August 13 (Barth et al. 
2003). The observations were made 4.65 h after the 
HETE-2 detection (Villasehor et al. 2002) when the vi- 
sual magnitude was ~ 19. The observed spectral range 
is A = 3220 — 9150 A, and the resolution is ~ 1 A per 
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Fig. 1.— Optical spectrum (AA = 1420 - 4080 A rest frame) of GRB 020813 (2 = 1.255; Barth et al. 2003) after correction for Galactic 
extinction with Eb-v = 0.111. The straight lino is the power-law best-fit F^ = i/-0.9i8±0.00l ^ iQia.SOSiO.Oll ^ 



Table 1 
Parameters for the absorption system at z 



1.255 



Ion 



log TV 

-2 



[cm-2] (km s-i) [X/Zn] 



Znll 


13.54 ±0.06 


62'' 




Sill 


16.29 ±0.04 


56 ±3 


-0.15 ±0.07 


Sill* 


14.32 ±0.11 


56^^ 




Mnll 


13.62 ±0.03 


81 ±27 


-0.66 ±0.07 


CrII 


13.95 ±0.03 


62*^ 


-0.62 ±0.07 


Fell 


15.48 ± 0.04 


62 ±1 


-0.91 ±0.07 


Nill 


14.20 ± 0.03 


43 ±5 


-0.94 ±0.07 


Till 


< 12.9^* 




< -1.0 


Mgll 


~ 15.8 


62'' 


0.7 


Mgl 


13.26 ±0.01 


110 ±6 




AlII 


~ 14.9*^ 


62'' 


0.4 


CI 


< 13. 7=* 






Call 


13.72 ±0.04 


. . .f 





''Effective Doppler parameter. 

''Doppler parameter is assumed to be the same as for Fe II. 

'^Doppler parameter is assumed to be the same as for Si II. 

4 (T upper limit, obtained assuming 6 = 62 km s~ . 
'^Column density estimated using the curve of growth. 
'^Ca II shows three absorbing components. 



pixel, or FWHM = 94-32 km s^^ from the blue to the 
red end. The signal-to-noise ratio in the ~ 3 hour total 
exposure is very good, S/N = 50 — 100, making this op- 
tical spectrum one of the best in terms of resolution and 
S/N ever obtained for a GRB afterglow. The spectrum 
of GRB 020813, after correction for Galactic extinction 
{Eb-v = 0.111; Schlegel et al. 1998), is shown in Fig- 
ure 1. Barth et al. (2003) reported the equivalent widths 
of detected absorption lines, and identified two absorp- 
tion systems, at z = 1.255 and z = 1.224 (see their 
Table 2). The former is very likely the redshift of the 
GRB. The lower-redshift system intersects the GRB line 



of sight, at a comoving distance^ of about 13 h~^ Mpc 
from the GRB. We fitted Voigt profiles to the absorp- 
tion lines to determine column densities for both absorp- 
tion systems, and cross-checked the results by using the 
curve-of-growth analysis (COG; Spitzer 1978) or the ap- 
parent optical depth method (Savage & Sembach 1991). 
As the intervening system at z = 1.224 has low metal 
column densities [^ 60 times lower than the z = 1.255 
system), we neglect its dust extinction ( ^ 0.01 mag in 
the optical). 

The absorption lines associated with the circumburst 
medium of GRB 020813 are numerous, and their equiv- 
alent widths Wr in many cases deviate little from the 
finear part of the COG (rest-frame Wr < 0.4 A), al- 
lowing a good determination of column densities. The 
main uncertainty is the total number of absorbing com- 
ponents, which is not known because of the finite reso- 
lution and the intrinsic blending of lines. Although the 
multiple absorbing components appear as a single line, 
the single-component approximation still gives robust re- 
sults (Jenkins 1986; Savage & Sembach 1991). How- 
ever, in this case, the width of the lines is an effective 
Doppler parameter and does not provide useful informa- 
tion on the temperature of the gas. Parameters of all 
lines relevant for this work (and 1 a errors) and abun- 
dances relative to zinc are given in Table 1. We adopt the 
common assumption, valid for neutral gas (A'h ~ -^m), 
that the abundances of ions with ionization potentials 
above 13.6 eV, are very close to the total abundance of 
the corresponding element: Nzn ~ -^Znii, -^Fe ~ -^FcII, 
Nq^ k, A'^crii, etc. We discuss possible contamination by 
high-ionization species (C IV and Al III) and ionization 
corrections in i l2.2l 

2.1. Low-ionization Absorption Lines 

Figure 2 shows, in velocity space, the best-fit Voigt 
profiles for the Fe II absorption in the GRB-DLA and 
in the intervening system at z = 1.224, which marginally 

^ We adopt cosmological parameters h = Ho/WO = 0.7, Cjv/ = 
0.3, and Qa = 0.7. 
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Fig. 2.— Absorption lines of Fe II in the GRB-DLA at z = 1.255 
(at km s~^) and of the intervening system at 2 = 1.224 (at —4200 
km s~^) and their best-fit absorption profiles {smooth lines). Nine 
Fe II absorption lines are detected in the GRB-DLA, allowing a 
good determination of the column density (log A'^pdi = 15.48 it 
0.04). In the top panel, we also show the best- fit Mn II triplet 
absorption (logAfMnll = 13.62 ±0.03). 



contaminates the GRB-DLA Fe II absorption. In the top 
and bottom panels, we also indicate the GRB-DLA Mn II 
triplet and the very weak Fe II A1611 line, respectively. 
Since many Fe II lines have large equivalent widths, we 
estimate the uncertainty due to saturation by consider- 
ing the weakest transitions, Fe II AA2249, 2260 [Wr < 0.3 
A). Using the COG, we derive logiVpcn = 15.5 ±0.1 and 
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Fig. 3. — Absorption lines of Zn II and Cr II. The smooth 
line is the best-fit absorption profile with logA^znii = 13.53 and 
logA^Crll = 13.95. The Doppler parameter is assumed to be the 
same as in Fe II (b = 62 km s~^). If the Doppler parameter 
is changed by ±15 km s~^, the best-fit Cr II and Zn II column 
densities change by ±0.01 and ±0.02 dex, respectively. The con- 
tamination of Zn II A2026 by Mg I A2026 is taken into account, and 
has been derived using the Mg I A2852 line fit (log A/^Mgl = 13.26). 
The Cr II A2026 absorption (not marked) is negligible. 



6 > 35 km s ^. This column density is consistent with 
(although less certain than) the one derived by profile 
fitting (logiVpoii = 15.48 ± 0.04 and 6 = 62 ± 1). Us- 
ing the apparent optical depth method, which gives reli- 
able results when the line profile is not known, we obtain 
logiVFeii — 15.47, in agreement with the profile-fitting 
result. 

In Figure 3 we show the profile fits for Zn II, Cr II, 
and Mg I. The Mg I A2026 line marginally contami- 
nates Zn II A2026 and is constrained by Mg I A2852 at 
A = 6433.4 A. The Doppler parameter b for Zn II and 
Cr II is very uncertain and thus has been assumed to be 
the same as for Fe II. Changing b by ±15 km s~^ does not 
have a significant effect on the column densities. In gen- 
eral, the column densities of lines with small equivalent 
widths depend only weakly on the Doppler parameter. 
Figure 4 shows, in velocity space, the best-fit Si II pro- 
files. Although the Si II A1526 line is contaminated by 
the z = 1.224 C IV AA1548, 1550 doublet, the fit is weU 
constrained by the blue side of the line and by the Si II 
A1808 absorption. 

In Figure 5, we show 4 weak {Wr < 0.2 A) Ni II lines, 
not identified by Barth ct al. (2003), together with the 
best-fit profile. In the observed wavelength range, there 
are also three Ti II transitions. Two are not detected, 
while a third, the stronger Ti II A3384 line, falls in the 
wavelength range of the telluric absorption at ^ 7600 A. 
However, we can still infer an interesting 4cr upper limit 
on the column density. 

The column densities of Mg II and Al II are very un- 
certain because these lines are heavily saturated. The 
presence of relatively strong Mg I and Ca II absorption 
(ionization potential below 8 eV) is another indication 
that the gas is mainly neutral. The non-detection of the 
relatively strong transition C I A1656 provides an upper 
limit for the column density. 

2.2. Higher Ionization Lines and the Ionization 
Correction 
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Fig. 4.— Absorption lines of Si II (logAfsiii = 16.29). The 
smooth Hne in the bottom panel also includes the best-fit absorp- 
tion profiles for the fine structure Si II* A1533 absorption and the 
C IV doublet of the intervening system at ^ = 1.224. We obtain 
a good fit to the Si II column density thanks to the isolated Si II 
A1808 absorption and the uncontaminated blue side of the Si II 
A1526 absorption {dashed line). 



The spectrum of GRB 020813 shows the higher ion- 
ization lines of C IV and Al III. We use these lines to 
demonstrate that the ionized gas is a small fraction of 
the total and that the ionization correction can be ne- 
glected. Al III and Al II absorptions were used to study 
the ionization state of the interstellar medium in QSO- 
DLAs (Howk & Sembach 1999). 

The C IV column density is very likely a small fraction 
of the total C column density in our GRB circumburst 
gas, because the product of the wavelength and the oscil- 
lator strength A/a for C IV A1551 is twice that for Fe II 
A2374, the cosmic abundance of C is 10 times higher than 
that of Fe, and C is usually much less depleted than Fe 
in dust ( < 10 times). Therefore, if C IV gas were a 
large fraction of the total C gas, we would expect the 
equivalent width of C IV A1550 to be much higher than 
that of Fe II A2374. However, the observed equivalent 
widths are 1.31 and 1.42 A for C IV A1548 and Fe II 
A2374, respectively (Earth et al. 2003), showing that the 
dominant ionization state of C is not C IV. Indeed we es- 
timate that the C IV column density is log A'civ — 13.5, 
compared with an expected total C column density (in 
gas form) of logiVc w 17. 

For Al, the ionization correction is less certain. The 
Al III doublet gives a column density of log A'^aiiii ~ 14.0 
(if its b value is the same as that of Fe II). Al II shows 
only the very strong line at A = 1670 A, for which we 
obtain logTVAiii ~ 14.9 (assuming b as in Fe II), thus 
Al III can only be a small correction of the total Al. If 
b is left as a free parameter, the Al III column density is 
even smaller, although not by much. We also note once 
more that most likely some Al III is located in a different, 
more ionized region in the host galaxy. 
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Fig. 5.— Absorption lines of Ni II (logAfNiii = 14.20). These 
weak lines have equivalent widths less than 0.2 A. If the Doppler 
parameter is assumed to be the same as for Fe II, the Ni II column 
density is 0.03 dex lower. 



2.3. The Si 11* Transition 

The spectrum of GRB 020813 also shows the ground 
fine-structure Si II* A1533 transition, from which we de- 
rive 



^Sill' ^ -[^Q-1.97±0.12 

iVsiii 



(1) 



The Si II* absorption was detected by Barth et al. (2003), 
but not identified. The iVsiii* /-^sm ratio can in princi- 
ple be used to determine the physical conditions of the 
gas (Bahcall & Wolf 1968). Theoretical predictions as a 
function of the gas temperature, gas density, and electron 
density have been made and compared with observations 
in QSO spectra (Srianand & Petitjcan 2000; Srianand & 
Petitjean 2001; Silva & Viegas 2002). The Nsiw/Nsm 
ratio we find in this GRB-DLA is ^ 100 times lower than 
the value measured in a QSO highly ionized region (Sri- 
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Table 2 
Heavy-element column densities 



logiVf 



Element log Af' 



obs 



WH^ 



WDH'" 



WD'^ 



CD^' 



Zn 

Si 
Mn 
Cr 
Fe 

Ni 



''Warm halo depletion pattern. 
''Warm disk+halo depletion pattern. 
•^Warm disk depletion pattern. 
'^Cool disk depletion pattern. 



13.54 ±0.06 13.53 ±0.06 

16.29 ±0.04 16.43 ±0.17 

13.62 ±0.03 14.27 ±0.03 

13.95 ±0.03 14.56 ±0.06 

15.48 ±0.04 16.38 ±0.19 



13.56 ±0.06 13.72 ±0.16 14.23 ± 0.09 

16.46 ±0.08 16.62 ±0.08 17.13 ±0.29 

14.30 ±0.08 14.46 ±0.05 14.97 ±0.17 

14.59 ±0.09 14.75 ±0.12 15.26 ±0.18 

16.41 ±0.12 16.57 ±0.10 17.08 ±0.12 



14.20 ±0.03 15.13 ±0.06 15.16 ±0.04 15.32 ± 0.03 15.83 ± 0.09 



7.5 



warm halo, D=1.05, xl =12.0 
warm disk+halo, D = 0.97, xl =67- 
warm disk, D = 0.96, xt =8-2 
cool disk, D=0.97, xl =35 
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Fig. 6. — Heavy-element column densities relative to solar abun- 
dances in the GRB 020813 circumburst medium. The dots are the 
measured values. The broken lines are the best-fit expectations 
if the Galactic dust depletion patterns are used as models. As a 
reference in the models, S is not locked into dust grains in all pat- 
terns, while Zn is largely depleted only in the cool disk clouds. The 
upper limit on Ti is shown but not used in the fit. The expected Ti 
values are from Savage & Sembach (1996) for the cool disk model, 
and are extrapolated for the other depletion models assuming that 
Ti scales as Ni (K. Sembach 2004, private communication). The 
measured Mg column density is shown {open circle) but not used 
in the best fit, because it is very uncertain. The dust-to-metals 
ratio, compared with the Galactic value D, and the best-fit xi a^r^ 
also reported. 



anand & Petitjean 2001). Although it is not possible to 
derive the density of the gas without knowing its temper- 
ature, the difference between the high-ionization gas in 
this QSO and GRB 020813 suggests a different physical 
nature. 

The Si II* transition was detected for the first time in a 
GRB afterglow by Vreeswijk et al. (2004), who reported 
a low value of the column density ratio, -/Vsiii/A^siii* ~ 
10"^'''. This GRB-DLA shows a very high H I column 
density (logA^ni = 21.9). Its gas is probably neutral, 
not only for the large H I column density, but also be- 
cause the Al III column density (estimated by us using 



the reported doublet equivalent widths) is more than 10 
times smaller than the expected total column of Al. We 
conclude that a low value of the A^siii*/-^Siii ratio (also 
found in our GRB-DLA) does not imply that the gas is 
ionized. The presence of Si II* is an important result 
that deserves further attention to understand better the 
physical conditions of the circumburst medium. 

3. DUST-DEPLETION PATTERN 

The absorption lines detected in the GRB spectrum 
give the column densities of heavy elements in the gas 
phase only. Thus, as different elements are locked in 
dust grains by different amounts, the column densities 
measured from absorption lines carry the signature of 
the chemical composition of the dust, i.e. the depletion 
pattern. 

We have accurate column density measurements of 
six heavy elements: zinc, chromium, silicon, manganese, 
iron and nickel, plus a stringent upper limit for titanium 
(Table 1). To derive the dust-depletion pattern, we com- 
pare these column densities with several dust-depletion 
patterns in the Milky Way (Savage & Sembach 1996): 
warm halo (WH), warm disk+halo (WDH), warm disk 
(WD), and cool disk (CD). 

We consider each of the four Milky Way (MW) deple- 
tion patterns as a model and denote by N^'^^ the ex- 
pected column density of clement i in such a model to be 
compared with the corresponding observed column den- 
sity N°'"' in the GRB spectrum. In our approach, N^^^ 
is also a function of the dust-to-metals ratio relative to 
the Galactic values {D = d/da)- 



N. 



exp 



N^'"[l + D{x,-l) 



(2) 



Here Xi is the Galactic fraction of element i in the gas 
phase and iV/°* is the total (gas and dust phase) col- 
umn density of the same element in the GRB-DLA. To 
express iVj*°', we define jji as the relative heavy-element 
abundance 

iV*°* = y, Ar™<=* , (3) 

where iV™*^* is the total column density of metals in 
the GRB-DLA. For convenience, we express yi in terms 
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Table 3 
Dust parameters 



Parameters 



WH^" 



WDHb 



WD<^ 



CDd 



D 



-Q,21J 



1.05 
-0.12 ±0.06 



0.97 0.96 0.97 

-0.09 ±0.06 0.07 ±0.08 0.58 ± 0.09 



0.40 



+0.06 



0.40 



+0.06 



0.56 



+0.11 
0.09 



1.84 



+0.43 



''Warm halo depiction pattern. 
''Warm disk+halo depletion pattern. 
'^Warm disk depletion pattern. 
''Cool disk depletion pattern. 



of solar abundances (Gravesse & Sauval 1998): i/i — 
c{Ni/N-ii)Q, where c is a constant. To obtain the best 
depletion model, we minimize the reduced x^ 



xl{D,N"^'^') 



1 



E 



logN""'' -log NiP{D,N' 



a{logN.^ 



obs\ 



(4) 

where ti = 6 is the number of elements in the fit and 
tT(logiV°'"*) is the error on the observed column density. 
We have defined the x" using logarithmic quantities, be- 
cause the observational errors have an approximate log- 
normal distribution. Using the observed column den- 
sities, we find four solutions (one for each of the four 
depletion patterns) with four different D and A^™''* val- 
ues. In Figure 6, we plot the column densities relative to 
solar abundances (both on a log scale) for the observed 
elements, together with the four best-fit models. The 
WH, WDH, and WD depletion patterns give very sim- 
ilar xt' whereas the CD pattern gives a worse fit. The 
results are similar if the errors on all the column den- 
sities are assumed to be equal. In Table 2, we report 
the total (gas+dust) column densities N-°* expected if 
the different depletion patterns in the MW are taken as 
models. The errors listed in Table 2 are the combination 
of 1 CT measurement errors and the deviations between 
the expectations given by the depletion models and the 
measured values. In Table 3 we present the overall dust 
parameters. 

4. EXTINCTION FROM DEPLETION PATTERN 

The extinction is proportional to the column of dust 
along the line of sight. In the solar neighborhood, a gas 
cloud with 10^"'^ cm~^ H I column density has Ay = 
0.5 visual extinction (Bohlin et al. 1978). The value 
of Ay is about 1/4 lower in a Large Magellanic Cloud 
(LMC) cloud with similar H I column density (Prevot 
et al. 1984); however, this is compensated by the lower 
metallicity (1/4 solar). Therefore, the visual extinctions 
for two clouds with similar columns of metals, one in the 
LMC and another in the solar neighborhood, are similar. 
The dust depletion in our GRB-DLA (where [Fe/Zn] = 
-0.91 ± 0.07) is closer to that in the LMC than that in 
the SMC ([Fe/Zn] = -1.12i^-^^ in the former, [Fe/Zn] 



= -0.57Io^? in the latter; Welty et al. 1997). 

If the dust in the GRB 020813 circumburst medium 
were like that in the MW or LMC, we would derive a 
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Fig. 7. — Abundances of Fe, Cr, Mn, and Si relative to Zn, 
vs. Zn II column density in QSO-DLAs {open squares) and 
in GRB-DLAs {filled dots). GRB-DLAs are: (a) GRB 990123, 
(fe) GRB 000926, (c) GRB 010222 (Savaglio et al. 2003), {d) 
GRB 020813 (this work), and (e) GRB 030323 (Vreeswijk et al. 
2004). For GRB 030323, the Zn II column density is not mea- 
sured, and we used S II instead, assuming [S/Zn] ~ 0. The mean 
errors for [X/Zn] in QSO-DLAs are ~ 0.08 dex. Dotted, short- 
dashed, long-dashed, and dot-dashed horizontal lines are the warm 
halo, warm disk+halo, warm disk, and cool disk values, respec- 
tively, reported by Savage & Sembach (1996). 
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Fig. 8. — Extinction curve index 7 and visual extinction Ay 
for constant confidence levels (CL), obtained by comparing tfie 
reddened intrinsic GRB emission with the observed GRB spectrum. 
The shaded vertical region indicates the Ay values not allowed by 
the depletion analysis (95% CL). 



visual extinction related to the column density of dust 
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Here A^^Hfi* and N^f^ 



are the dust and metal columns, 



respectively, for a A'hi = 10 cm^ gas cloud with solar 
metallicity, and D and iV""^*, as before, are the dust-to- 
metals ratio and the total metal column in our GRB- 
DLA. The A^™'=77V™|*i and Ay values (and relative 1 a 
errors) are listed in Table 3. The inferred Ay values are 
all higher than 0.3 (95% confidence level), with a best 
fit, obtained for the WDH depletion pattern (Fig. 6) of 
Ay = 0.40 
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4.1. QSO-DLAs and GRB-DLAs 

The QSO-DLAs are usually considered to be represen- 
tative of the ISM of the whole population of galaxies 
at high redshifts. The dust effects in QSO-DLAs are 
generally small: the visual extinction in 44 QSO-DLAs 
(with detected Zn II, 0.39 < z < 3.39) inferred using 
the same approach as above, is always Ay ^ 0.15. This 
difference from the extinction in the GRB 020813 cir- 
cumburst medium {Ay > 0.3) is indicated by the lower 
Zn II column densities and the larger [Fe/Zn] values typ- 
ically seen in QSO-DLAs (zinc and iron are among the 
least and most depleted elements, respectively). This is 
illustrated in the bottom panel of Figure 7. 

The difference between GRB-DLAs and QSO-DLAs 
may be the result of an observational bias: QSO-DLAs 
are generally found in the foreground of optically se- 
lected QSOs, and dusty QSO-DLAs are hard to find in 
such samples (Fall & Pei 1993). GRBs are (at least for 
some time) much brighter than QSOs and are detected 
by gamma-rays and so might be much less affected by 
this bias. Moreover, we know that GRBs are located in 
star-forming regions, where the densities are presumably 
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Fig. 9. — fntrinsic and reddened GRB spectrum {upper and 
lower dotted lines) as compared with the observed spectrum {solid 
line), for different extinction laws, fn all cases, larger Ay values 
are not allowed (95% CL). Note that the intrinsic GRB emission is 
obtained by the best fit and is not assumed, fn the top panel, the 
extinction curve is a power law where 7 = 0.845 is the maximum 
allowed value for Ay = 0.40. In the bottom three panels, MW, 
LMC, and SMC are assumed. 



higher than average, and large column densities of dust 
are not surprising. 

The elements Cr, Mn, and Si (in addition to Fe) also 
indicate relatively high dust depletion in GRB 020813 
(Fig. 7). Similar high depletions were found in the cir- 
cumburst medium of three other GRBs (Savaglio et al. 
2003). In another GRB-DLA (in GRB 030323), the 
[Fe/S] is consistent with zero, indicating little or no dust 
depletion (Vreeswijk et al. 2004). Since Zn II is not mea- 
sured in this object, while S II is, we assume [Zn/S] ~ 
and report this point in Figure 7. 

5. EXTINCTION FROM CONTINUUM SPECTRUM 

The presence of dust in the GRB circumburst medium 
can affect the observed GRB continuum spectrum. The 
GRB extinction law has never been determined by ob- 
servation and is not known a priori. From the lack of 
apparent reddening in many GRB spectra, it has been 
concluded in those cases that dust is absent or that the 
extinction law is a weak function of the wavelength (gray 
extinction; Stratta et al. 2004). Perna et al. (2003) have 
shown, using numerical simulations, that the GRB X- 
ray/UV radiation field can significantly flatten the dust 
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extinction curve. 

We have already derived the heavy-element depletion 
pattern in GRB 020813, from which we inferred the vi- 
sual extinction to be relatively high, Ay > 0.3, assuming 
the same rate of visual extinction as in the MW or LMC 
per unit of metal column density. However, the observed 
spectrum shows very little deviation from a perfect power 
law and no sign of the 2200 A extinction feature (Fig. 1), 
which allows us to place strong constraints on the ex- 
tinction curve of the intervening dust. We use a fairly 
general approach, namely, we assume that the intrinsic 
GRB emission is a power law 



n-^i^y 



(6) 



where Fy and a are not fixed a priori, and we apply a 
dust extinction law of the form 



. /5500V 



(7) 



We vary a, Fy, 7 and Ay^ and compare the correspond- 
ing model spectrum with the observed GRB 020813 spec- 
trum using x^ minimization. Because we are mainly 
interested in the confidence levels (CL) of only two of 
these four parameters [Ay and 7), we project the four- 
dimensional confidence regions in the Ay and 7 space. 

In Figure 8 we show 7 versus Ay for constant CL [solid 
lines). For Ay = 0.40 (as inferred from the dust deple- 
tion), any value of 7 between 0.85 and is allowed (95% 
CL), and this would give an intrinsic GRB spectral slope 
a < 1.72. The shaded vertical region is the one excluded 
by our depletion-pattern analysis. 

We also calculated the maximum Ay values (95% CL) 
obtained with MW, LMC, and SMC extinction curves. 
For the MW and LMC extinctions, these give Ay < 0.08 
and < 0.19, respectively, resulting from the absence of 
the 2200 A extinction feature in the observed GRB spec- 
trum. For the SMC extinction, which has a very weak 
2200 A bump, the constraint is Ay < 0.18. 

In Figure 9, we show the intrinsic and reddened GRB 
model spectrum as compared to the observed spectrum. 
For the MW and MC extinction curves, this is done for 
the maximum allowed Ay value (bottom three panels), 
while for the power-law extinction [top panel), Ay = 0.40 
is assumed. Since the upper limits on the visual extinc- 
tions inferred from MW and MC reddening {Ay < 0.2) 
arc smaller than the lower limits inferred from the de- 
pletions {Ay > 0.3), we have reached a contradiction 
and are forced to conclude that the dust in the circum- 
burst medium of GRB 020813 cannot be like that in the 
MW and MC in all respects (both chemical and optical 
properties). We note that although variations in the ex- 
tinction law in different OB association stars have been 
found in the MW (Fitzpatrick & Massa 1988), the pres- 
ence of the 2200 A bump is still very clear and would not 
be consistent with the observed GRB spectrum. 



6. SUMMARY AND CONCLUSIONS 

In this work, wc studied, for the first time simultane- 
ously, the dust depletion pattern and the extinction law 
in the circumburst medium of GRB 020813. Very little 
is known about the dust properties in the circumburst 
medium of any other GRB, mainly because the emission 
flux declines so rapidly that optical spectra are very hard 
to obtain. 

An excellent optical spectrum of the GRB 020813 af- 
terglow {z ~ 1.255) was obtained with the Keck I LRIS 
'--^ 5 hr after the burst by Barth et al. (2003), when the 
GRB was still rather bright (19th mag in the optical). 
From the numerous and strong absorption lines detected 
in the spectrum, we measured the column density of six 
heavy elements (Zn, Si, Cr, Mn, Fe, and Ni). The relative 
abundances resemble the Galactic dust-depletion pat- 
terns, from which we infer a visual extinction Ay ~ 0.40, 
and Ay > 0.3 (95% confidence). 

The relatively high visual extinction in GRB 020813 
is hard to reconcile with the observed UV GRB contin- 
uum spectrum, which deviates very little from a per- 
fect power law, suggesting low reddening. However, 
this can be explained if the extinction curve is only a 
weak function of the wavelength. We constrained the 
shape of the extinction curve, assuming the power-law 
form Ax = Ay(5500/A)T. For Ay = 0.40, we obtain 
7 < 0.85. This high visual extinction rules out the pos- 
sibility of MW, LMC, or SMC reddening, for which only 
Ay < 0.08, 0.19, 0.18, respectively, are allowed. The MW 
and LMC extinctions are also constrained by the absence 
of the 2200 A extinction feature. 

A moderately high visual extinction is not surprising if 
GRBs originate in regions of intense star formation and 
high density (Hjorth et al. 2003a). However, under some 
conditions, the GRB itself could destroy some of the dust 
grains (Fruchter et al. 2001). 

Finally, heavy-element column densities in the circum- 
burst medium of this GRB are considerably higher than 
those typically observed in QSO damped Lya systems 
(about 10 times higher), which is likely the result of a 
selection bias in the latter. The comparison of the dust- 
depleted element Fe, with the non-depleted element Zn, 
in the gas phase, gives [Fe/Zn] — —0.91 ±0.07, whereas in 
44 QSO-DLAs {{z) = 2.0) this is [Fc/Zn] = -0.53 (0.28 
dex 1 a dispersion), suggesting a higher dust content in 
the circumburst medium of GRB 020813. 
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